Kurtosis and compressibility near the chiral phase transition 
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The properties of net quark number fluctuations in the vicinity of the QCD chiral phase transition 
are discussed in terms of an effective chiral model in the mean-field approximation. We focus on 
the ratio of the fourth- to second- order cumulants (kurtosis) and the compressibility of the system 
and discuss their dependence on the pion mass. It is shown that near the chiral phase transition, 
both observables are sensitive to the value of m,r. For physical m^, the kurtosis exhibits a peak 
whereas the inverse compressibility shows a dip at the pseudocritical temperature. These structures 
disappear for large m-n-. Our results, obtained in an effective model with two flavors, are qualitatively 
consistent with recent results of 2-\-l flavor lattice gauge theory. We also discuss the high- and low- 
temperature properties of these observables and the role of the coupling of the quark degrees of 
freedom to the Polyakov loop. 

PACS numbers: 



I. INTRODUCTION 

Understanding the critical properties of the QCD 
medium and its phase diagram is a central problem in 
strong interaction physics, which is being addressed in 
both theoretical and experimental studies. Lattice Gauge 
Theory (LGT), in particular, has provided substantial in- 
sight into the physics of the transition between the low 
temperature hadronic and high temperature quark-gluon 
plasma phase [H, i, i, i, H, 11, 0, S, H- Effective models 
also provide a useful qualitative description of the chiral 
phase transition [l^, [ll|, [H, [H, Here, the expecta- 
tion that the critical behavior of such models is governed 
by the same universality class as QCD [l^ is of cen- 
tral importance. Furthermore, recent experimental data 
from heavy-ion collisions at RHIC and SPS provide new 
insight into the properties and behavior of a deconfined 
medium created in dynamical systems [T6j . 

The phase diagram of a thermodynamic system may be 
mapped out by studying fluctuations and their response 
to changes of the thermodynamic parameters. In QCD 
the fluctuations of conserved charges are particularly per- 
tinent for exploring deconfinement and chiral symmetry 
restoration @, 0, M, [H, [H, US ■ 

Studies of QCD on the lattice at finite temperature and 
density show that the ratio of the quartic to quadratic 
fluctuations of the net quark number Ri 2 (kurtosis ^) is 
a valuable probe of deconfinement and chiral dynamics 
0,11, Hi- Furthermore, both model calculations [2ll. [2^ 



Strictly speaking, in statistics kurtosis is given by the 
fourth moment divided by the second moment squared, i.e., 
{(SN)^)/ {{SN)'^)^ . In this letter we employ the so called kurtosis 
excess, with an atypical normalization i?4^2 = ((^-'V)*)/((5Af)^) — 
3((i5A'^)^). It turns out that this ratio is useful for distinguishing 
the confined and deconfined phases 0| . In the interest of brevity, 
we sacrifice rigor and refer to -R4,2 as kurtosis. 



and LGT studies 0, 0j @| indicate that the inverse com- 
pressibility i?K is a useful observable for identifying the 
position of the critical end point CEP in the QCD phase 
diagram. 

In the limits of high and low temperatures, i?4^2 de- 
pends on the net quark content of the relevant baryon 
number carrying degrees of freedom [3]. Consequently, 
one expects a strong reduction of i?4^2 at the confinement- 
deconfinement transition. Furthermore, the inverse com- 
pressibility Rk, vanishes at the CEP, due to the divergence 
of the net quark number fluctuations at a critical point, 
belonging to the Z{2) universality class The first 

LGT studies of these observables, which were performed 
in 2-fiavor QCD with a pion mass of ~ 770 MeV 0,01, 
confirmed the expected behavior of R4 2 above and be- 
low the deconfinement transition. However, the results 
for the inverse compressibility R^ shows no sign of the 
expected decrease at the crossover temperature Tpc] the 
LGT results for both i?4.2 and i?„ interpolate smoothly 
between the limiting values in the hadronic and quark 
gluon plasma phases. 

By contrast, results for i?4_2 and i?„ obtained recently 
in 2+1-flavor QCD with an almost physical pion mass, 
jtItt — 220 MeV, exhibit a qualitatively different temper- 
ature dependence @. In these calculations, the kurtosis 
was found to develop a peak whereas the inverse com- 
pressibility acquired a dip at the pseudocritical temper- 
ature. These results indicate that, close to the phase 
transition, these ratios depend on the strength of the ex- 
plicit symmetry breaking term, i.e. on the value of the 
pion mass. 

In this letter, we use universality arguments, relating 
QCD with chiral effective models, to explore the proper- 
ties of the kurtosis and the inverse compressibility in hot 
and dense matter. We focus on the dependence of i?4^2 
and Rk on the pion mass near the chiral transition. As 
an effective model of QCD we use the quark-meson La- 
grangian coupled to the Polyakov loop in the mean field 
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approximation [l4| . 

We show, that the model results for R4 2 and 
are consistent with the LGT results and find that the 
peak/dip structures are indeed strongly dependent on the 
value of the pion mass. Given this, we conclude that the 
non-monotonic behavior of these ratios found in LGT, for 
an almost physical pion mass, is due to the 0(4) chiral 
dynamics of 2-flavor QCD. Since the explicit symme- 
try breaking in the strange sector is large, one expects 
the strange quarks to decouple from the critical fluctu- 
ations and consequently that the critical properties, for 
physical quark masses, are dominated by the up-down 
sector. This observation supports the notion that the 
scaling properties of 2-|-l-flavor QCD are also governed 
by the 0(4) universality class. 



II. THE CHIRAL MODEL 

In order to explore the influence of chiral dynamics on 
the above observables we use the effective Lagrangian of 
quark-meson model coupled to the Polyakov loop (PQM 
model) [3 



The gluonic contribution to the thermodynamic po- 
tential is parametrized as an effective potential for the 
Polyakov loop and its conjugate 



where 
62 (T) 



2 £*£-^{e^ + r^) + ^{e*e)\ (2.5) 



02 



as 



.(2.6) 



The temperature independent coefficients, at and bi in 
Eqs. (|2.5p and (|2.6p have the following values: qq = 6.75, 
ai = -1.95, 02 = 2.625, = -7.44, 63 = 0.75 and 
64 = 7.5. They are chosen so that the LGT results for 
the equation of state of a pure SU{3) gauge theory are 
reproduced [lol |. Furthermore, for the parameter To we 
choose To = 208 MeV, the value favored in ref. for 
the two-flavor sector. 

The quark/antiquark contribution to the thermody- 
namic potential is given by 



-P,,{T,y.)^-2NfT 



(2^ 



(2.7) 



£ = q{^P-g{a + ^J,rn))q+^-{^,ar + l{^,nr {ln[l + 3(^ + re-(^-'')/^)e-(^-^)/^4.-3(^-'')/^] + 

(2.1) ln[l + 3(r + £e^(^''+^)/^)e-(^-+'^)/^^-3(^-+'')/^] } , 



with U{i, £*) being an effective potential of the gluon field 
expressed in terms of the traced Polyakov loop £ and its 
conjugate £*. The 0(4) representation of the meson fields 
is = (cr, 7?) and the corresponding S'C/(2)L(3'5'[/(2)fl chi- 
ral representation is given hy a + if ■ nj^. Consequently, 
there are 7V| = 4 mesonic degrees of freedom coupled 
to Nf = 2 flavors of constituent quarks q. The effect 
of gluons are accounted for by introducing an effective 
field, which couples to the quarks by means of a covari- 
ant derivative 



where /i is the quark chemical potential and Ep 



(2.2) 



Here the spatial components of the gauge field are chosen 
such that their average vanish, i.e. = 5^o^o- 

In the mean field approximation the thermodynamics 
of the model is described by the thermodynamic potential 



p'^ + is the quasi-particle energy with rriq = ga 

being the constituent quark mass. 

In the mean field approximation the three independent 
fields cr, I and I* are approximated by their thermal av- 
erage values, which are determined by extremizing the 
thermodynamic potential. 



III. KURTOSIS AND COMPRESSIBILITY 
NEAR THE CHIRAL TRANSITION 

In the following we focus on the quark sector of the 
thermodynamic potential. Following the procedure in 
LGT studies we use the Taylor expansion in /i/T to 
compute various observables. In particular the thermo- 
dynamic pressure is expressed as 



n = U{£J*) + U{a)+ngg{£J*,a), (2.3) 



with the gluonic, mesonic and quark/antiquark contribu- 
tions, respectively. 

The mesonic part of the potential is of the form 

U{a, Tf)^^ {a^ +TT^-v^Y - ha, (2.4) 

where the last term describes the interaction of the scalar 
field a with an external field h. This term explicitly 
breaks the chiral symmetry of the Lagrangian and gives 
the pion a finite mass = h/a. 



rp4 



with the expansion coefficients 

9"(p(T,/i)/T4) 



c„(T) 



(3.1) 



(3.2) 



The Taylor coefficients are generalized susceptibilities 
corresponding to moments of the net quark number 
QjHjHI- In particular, the first two non- vanishing deriva- 
tives, C2 and C4 are the second and fourth order cumulants 
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FIG. 1: The kurtosis i?4,2 for various values of the pion mass in the PQM (left panel) and in the chiral quark-meson (QM) 
model (right panel). The temperature is normalized to the pseudocritical temperature Tc for a physical value of the pion mass. 



which are related to fluctuations of net quark number, 
5Nq ^ Nq- (Ng), as follows: 



Thus, in the limit of high temperatures, i?4,2 and are 
in both models given by 



C2 



C4 



^ =< iSN,^' 



< Nq 



< i6Nq)^ > -3 < (SNq)' >2 



(3.3) 



with Xq being the quark number susceptibility. We focus 
on the ratio 



R - ^4 

C2 



(3.4) 



which we define as the kurtosis of the quark number fluc- 
tuations (see footnote above). We consider the dimen- 
sionless ratio of the inverse compressibility and Gibb's 
free energy density ^Ug, 



Rk. — 



(3.5) 



where the inverse isothermal compressibility of the sys- 
tem is given by 



■^4,2 



T»To 



R> 



«IT>>To 



3(m/T)2 



(3.8) 



where i?4^2 is evaluated at ^ = 0, while for R^ we re- 
tain the fi dependence in anticipation of the application 
below. 

One of the main features of a chiral Lagrangian cou- 
pled to the Polyakov loop is the " statistical confinement" 
which implies, that at small T, the effective degrees of 
freedom are three quark states [13, [Hi . This is explicitly 
seen in Eq. (|2.7[) . where for small T «To the Polyakov 
loop / ~ 0, implying suppression of the one- and two- 
quark states in the partition sum. Consequently, in the 
Boltzmann approximation" the fermion contribution to 
the pressure is 



PqqiT,^l) 2Nf (im 



rp4 



277r2 



T 



3mg 



cosh|^. (3.9) 



-V 



dp 
dV 



,uq 

Xq 



(3.6) 



Before considering i?4_2 and i?„ near the chiral phase 
transition, we first focus on their behavior well above 
and below the transition region. For high temperatures, 
i.e. T » To, the ratio niq/T « 1 and the Polyakov 
loop / — > 1. Hence, the quark- antiquark contribution to 
the thermodynamic pressure (|2.7p converges to that of 
an ideal gas of quarks, where 



PqqiT,fi) 
2^4 



NfNc 



1 



127r2 



7f_ 
180 



(3.7) 



Thus, at low temperatures the quark/antiquark contri- 
bution to the thermodynamic pressure is that of a non- 
interacting gas composed of particles and antiparticles 
with mass M = 3mq and baryon number B = 1 and 
B — respectively. The effective degeneracy of these 
particles is 2Nf/27. The suppression by a factor 1/27 
is due to the rescaling of the momenta implied by the 
transformation 3\/p^ + -^/fc^ + (3m)2. 



^ The Boltzmann approximation is valid if Sniq/T > 1. This con- 
dition is clearly satisfied for T < To 
MeV. 
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An important feature of the pressure (|3.9p is the fac- 
torization of the MjT and /i/T dependence, 



(3.10) 



A similar factorization occurs in the hadron resonance 
gas model HRG, although with a more complicated func- 
tion replacing f{M/T). In the HRG, the corresponding 
function involves a sum over all baryon states 0,111, H^l ■ 
We note, however, that the Polyakov-loop models fail to 
reproduce the pressure of a resonance gas, even when the 
HRG is restricted to a gas of nucleons and A resonances 
(more generally the lowest three-quark multiplets) and 
their masses are approximated by Mn ~ Ma ~ Sm^. 
The origin of this failure is twofold: on the one hand the 
suppression factor mentioned above and, on the other 
hand, an incorrect dependence on Nf. The latter is due 
to the fact that only states where all three quarks have 
the same flavor are counted in p.9p . 

Since the HRG seems to provide a satisfactory de- 
scription of the QCD equation of state at low tempera- 
tures 0, m, [131 , the effective model considered here can, 
strictly speaking, be valid only near the phase transition, 
where the dynamics is controlled by the chiral symmetry. 
However, the model may also be useful outside the criti- 
cal region, for computing observables that are not sensi- 
tive to details of the mass spectrum. The kurtosis p.4p 
and the dimensionless inverse compressibility (|3.5p fulfil 
this criterion at least at low temperatures. Indeed, from 
Eqs. (P^ . ([X5|l and ([XTI])) it follows that the function 
/(A//T), which depends on the mass spectrum, cancels 
in the ratios, leading to 



i?4 



2It<<To=9' -^'ilT<<Tr, - — tanh^^. (3.11) 



3/i 



T 



For reference we also give the dimensionless inverse com- 
pressibility in the approximation, where Uq and Xq ^-rc 
computed to next-to- leading order in ^/T 



_ 2 + 3(m/T)^ 



(3.12) 



The kurtosis (|3.4p is strongly dependent on the net 
quark content of the baryon number carrying effective 
degrees of freedom. For instance, at low temperatures, 
Ra,2 = (35)^. Consequently, if the low-temperature 
phase is confining, B = 1 and i?4^2 = 9, as found also 
in the hadron resonance gas model 0, 0] ■ On the other 
hand, in the pure Quark-Meson (QM) model, i.e. neglect- 
ing the effects of the Polyakov loop, the relevant fermionic 
degrees of freedom are single quarks, at all temperatures. 
Therefore, the corresponding value of the kurtosis in the 
low temperature limit is Ra,2 = 1- Clearly, in an ef- 
fective chiral model with quark degrees of freedom, the 
Polyakov loop is essential for obtaining the correct low 



temperature behavior of the kurtosis ^ . 

We now turn to the temperature dependence of i?4 2 
and Rk- By inspection of Eqs. p.8p and (|3.1ip we con- 
clude that both the kurtosis and i?^ are, for fixed /i/T, 
temperature independent far above and below the transi- 
tion. However, the high- and low-temperature values dif- 
fer considerably, implying that the two ratios must vary 
near the phase transition. The major contribution to the 
change of the kurtosis i?4,2 between the low and high 
temperature limits is due to the change of the relevant 
baryon number carrying degrees of freedom, from three- 
quark states to single quarks. Thus, indeed the kurtosis 
is, as noted in 0|, an excellent probe of deconfinement. 
It is clear from this discussion that, far above and be- 
low the transition, i?4_2 and i?^ are independent of m-r^. 
However, as illustrated in Fig. [1] this is not the case in 
the phase transition region. 

In Fig. [1] we also show, that for T < Tpc the kurtosis in 
both the PQM and QM models converges fairly rapidly 
to its low-temperature limit, independently of the value 
of ttItt. The temperature (in)dependence of i?4.2 in both 
models is consistent with LGT findings in 2- and 2-1-1- 
flavor QCD 0, 0. However, as discussed above, only 
the PQM model correctly reproduces the low T limit, 

i?4,2|T«Tp, = 9- 

At high temperatures, the kurtosis in the QM model 
reaches the asymptotic value already at T ~ l.STpc, hi- 
dependently on the value of the pion mass, in full agree- 
ment with LGT results 0, Q. However, this is not the 
case for the PQM model, where even at T > 2Tpc, i?4,2 is 
still about 50% above the ideal gas value. The different 
high-temperature behavior of the PQM and QM models 
is clearly a consequence of the coupling of the quarks to 
the Polyakov loop rather than an effect of the mean field 
dynamics, as conjectured in Ref. [l^. The ideal gas value 
of i?4_2 is in the PQM model reached at high T if both 
rUq/T << 1 and ^ ~ 1. The first requirement is indeed 
satisfied already at T > (1.2 — 1.5)Tc, but the second is 
not valid even at T > 2Tpc- This is a direct consequence 
of the parametrization of the Polyakov loop effective po- 
tential in the PQM Lagrangian (|2.ip . The parameters 
of the effective potential were fixed by fitting SU (3) lat- 
tice results, where I approaches its asymptotic value very 
slowly. This discussion points to a clear limitation of the 
applicability of the PQM model to the interpretation of 
LGT results at high temperatures where e.g., the interac- 
tions of quarks with space- like Wilson loops, not included 
in the PQM model, could be relevant. 

Near the chiral transition, the kurtosis is strongly de- 
pendent on the pion mass. For a physical pion mass, i?4^2 
exhibits a peak in both models. With increasing the 



^ We stress again that Polyakov loop models can reproduce the 
correct value of Ra,2, although they do not yield a realistic de- 
scription of the pressure at low temperatures, because the de- 
pendence on details of the baryon spectrum cancels in the ratio 
C4/C2. 
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T/T, T/T^ 

c c 

FIG. 2: The inverse dimensionless compressibility computed in the PQM model at fixed n/T = 0.8 for various values of 
rrin (left panel) and at fixed = 138 MeV for different fi/T (right panel). The curve labelled CEP corresponds to the value 
of /i/T at the critical end point. The arrows in the right panel indicate the low-temperature values of ii^ at a given /i/T. The 
temperature is normalized as in the Fig. [1] 



peak height gradually decreases and its position is shifted 
to larger temperatures, reflecting an increase in Tpc- The 
dependence of the kurtosis on the pion mass is stronger 
in the QM than in the PQM model. A comparison with 
LGT results 0, Q show that the PQM model yields a 
much better overall description of the LGT data. For 
large 171-^ there is a smooth change of the kurtosis be- 
tween its limiting high- and low-tcmpcrature values, in 
agreement with 2-flavor QCD on the lattice 0, 01 ■ On 
the other hand, for physical m^r, i?4.2 increases at Tpc 
beyond the low-temperature value i?4_2 — 9. Both these 
results, together with the low temperature limit, are con- 
sistent with the LGT results in 2 and 2-1-1 flavor QCD, 
obtained for — 770 and TOtt ~ 220 McV, respectively. 
The only problem of the PQM model in this context is 
the high-tcmperature behavior of i?4.2, discussed above, 
while the QM model yields an incorrect low-temperature 
limit and a too strong dependence on the pion mass near 
the phase transition. Consequently, we restrict our dis- 
cussion below to the PQM model. 

As a consequence of 0(4) scaling, one expects a strong 
variation of the kurtosis with m^r near the chiral phase 
transition 

0. At /X = 0, the singular part of the net 
quark number fluctuations scales as Xg ~ t^~°' whereas 
that of the fourth order cumulant is of the form C4 ^ t~°', 
where t = \T—Tc\/Tc and a is the critical exponent of the 
specific heat. In the 0(4) universality class, a is small 
and negative, a c± —0.26. Consequently, the regular part 
dominates in the susceptibility Xqi whereas the singular 
part of C4 corresponds to a cusp. In the presence of a 
term that explicitly breaks the 0(4) symmetry, the cusp 
is smoothened and can, for small pion masses, be seen in 
i?4_2 as a peak at Tpc- For larger ttItt, C4 is dominated by 
the regular part, resulting in monotonic change of i?4.2 



in the transition region. These expectations, which are 
borne out in the PQM model, are consistent with LGT 
results. 

At finite chemical potential, the singularity is stronger; 
C4 diverges at the 0(4) line whereas the quark number 
susceptibility Xg develops a cusp. Furthermore, Xg di- 
verges at the critical end point (CEP), in accordance 
with Z{2) universality. Consequently, the inverse com- 
pressibility, introduced in the Eq. (j3.5p . is an interesting 
observable that can be used to verify the existence of the 
CEP, to identify its position, provided it exists and to 
establish the universality class of the chiral transition. 

In left panel of Fig. [2] we show the temperature de- 
pendence of the dimensionless inverse compressibility 
([XS]) computed in the PQM model for fixed fx/T = 0.8 
and for various values of . We follow the procedure in 
the LGT calculations of 0, H i and compute Uq and Xg 
keeping only the first two non-vanishing coefficients, C2 
and C4. The low T values of i?„ are consistent with the 
low-temperature limit (j3.12p . indicated by the arrows in 
the right panel. For large temperatures, T » Tpc, is 
independent of m^r and slowly converges to the ideal gas 
value p.Sp . Similarly as in the case of i?4,2, Rk deviates 
substantially from the Stefan-Boltzmann result, due to 
the gentle approach of the Polyakov loop to unity. 

Near the chiral transition, i?„ shows a rather strong de- 
pendence on the pion mass. For large rriTr, Rk is smoothly 
increasing between the low- and high-temperature values. 
On the other hand, for a physical value of the pion mass 
(771^ = 138 MeV), a dip develops in the vicinity of Tpc- As 
shown in the right panel of Fig. [2l the value of R^ at the 
dip drops with increasing /i/T. This is consistent with 
recent lattice results Q obtained in 2-1-1 fiavor QCD for 
TOtt = 220 MeV. The agreement of the two-flavor PQM 
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model results with those of 2+1 flavor LGT suggests that 
the lattice results reflect the critical behavior of 0(4) chi- 
ral dynamics. 

For the parameters used in our calculations, the CEP 
appears in the PQM model at /i/T ~ 1. Consequently, at 
this value of /i/T, Rk should vanish at the critical tem- 
perature. It is clear from the right panel of Fig.[2]that 
clearly reflects the 0(4) critical dynamics. However, the 
approximate i?^, obtained by a Taylor expansion of the 
numerator and the denominator, is not adequate to iden- 
tify the position of the CEP by a zero of Although 
the value of R^ at Tpc drops with increasing n/T, as the 
CEP is approached, it remains different from zero. This 
is clearly a consequence of the Taylor expansion, since the 
zero of Rk is due to the singularity of the quark number 
susceptibility, which is not reproduced by a polynomial 
of finite order. 



IV. CONCLUSIONS 

Based on the quark-meson chiral Lagrangian coupled 
to the Polyakov loop we have explored the properties of 
the kurtosis i?4.2 of the net quark number fluctuations 
and the dimensionless inverse compressibility We 
have focused on the behavior of these observables near 
chiral transition as functions of the pion mass. 

We have shown, that it is necessary to include the 



coupling of the quarks to the Polyakov loop, in order 
to reproduce the low temperature behavior of the kur- 
tosis expected in QCD. We have also demonstrated that 
near the pseudocritical temperature, both these observ- 
ables are very sensitive to the pion mass. For a physical 
value of TOtt , the kurtosis develops a peak and the inverse 
compressibility a dip at the pseudocritical temperature. 
These structures disappear for large mT^. Our results 
supports the notion that, for physical quark masses, the 
strange quark decouples from the critical fluctuations and 
the scaling properties are governed by 0(4), the expected 
universality class of two flavor QCD. We also found a 
suppression of i?^, as the critical end point CEP is ap- 
proached. However, i?„ remains finite at the CEP, if the 
Taylor expansion in /i/T of the thermodynamic pressure 
is implemented. 
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